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ABSTRACT 

The Milky Way spiral arm dynamics is studied using the birth place of open clusters. We sep¬ 
arately study the nearby spiral arms, and find evidence for multiple spiral sets. In particular, 
the Sagittarius-Carina arm appears to be a superposition of two sets. The first has a pattern 
speed of Op^carina.i = ± l.lstat km sec“^ kpc“\ while the second has Dp,carina ,2 = 

29.8ti',3sys i l-3siat km sec“^ kpc“^. The Perseus arm located at larger galactic radii exhibits only 
one clear set with Op^Perseus = ‘20-0t.\'^sysH-6stat km sec“^ kpc“^. The Perseus and “slower” Carina 
arms are most probably part of an m = 4 set. The Orion “armlet” appears to have a pattern speed 
of Op^orion = ± O.Sgtat km sec“^ kpc“^. In other words, both the “faster” Carina arm and 

the Orion arm are nearly corotating with the solar system. 

Subject headings: Galaxy: kinematics and dynamics, open clusters and associations: general 


1. INTRODUCTION 

There is currently no doubt that the Milky Way is 
a spiral armed galaxy. Yet, our edge-on view and dust 
obscuration through the galactic disk do not permit us to 
see the structure and pattern of the spiral arms directly. 
Thus, different claims for different spiral arm geometries 
and dynamics are abundant in the literature. Claims 
for a double or 4-armed structure are common, while 
numerous estimates for the pattern speeds can be found. 

A 4-armed spiral pattern does however ap- 
pear to be preferred by the observa tional data 
i|Georgelin fc GeorgeliiJ Il976t IValled 1199511. In par¬ 
ticular, velocity-longitude maps reveal quite clearly the 
existe nce of a 4-armed spiral pattern outs ide the solar 
circle l|B1itz et al\ 11 98, 'll: iDa.me et ~ai\ I 2 OOIU . Unfortu¬ 
nately, however, nothing equivocal can be stated at 
smaller galactic radii because of the inherent ambiguity 
in V — £ maps and the large scatter in the indirect tracers 
of spiral arms. 

Besides the number of arms m, in the spiral set, the spi¬ 
ral arms are also characterized by a pitch angle i, a pat¬ 
tern speed Hp and an amplitude. For a logarithmic arm, 
the pitch angle is fixed and given by cot i = \R{d(j)/dR )\, 
where R{(j)) is the parametric description of the ar m. By 
consid ering the various observations at the time, iVafled 
1 ( 19 ^ has shown that an m = 4 set, with a pitch angle 
of z = 12 ± 1°, is notably more consistent with the data. 

Notably less consensus, however, can be found for the 
value of the pattern speed of the spiral arms. Estimated 
values, based on different methods and data sets, typi¬ 
cally range bet ween 10 to 30 km sec“^ kpc“^ (for a re¬ 
cent review, see lShavi^l200,‘l() . but the distribution is not 
even—the ~ 20 results seem to divide into three different 
ranges of values. These are ~ 9 — 13, ~ 1 — 5 

and ~ (—1) — (—4) km sec“^ kpc“^. 

We can classify the methods according to the core 
assumptions upon which they are based. In particu- 

^ The maps reveal 3 arms separated by ~ 90°, a fourth arm is 
presumably obscured behind the center of the galaxy. 


lar, several methods implicitly assum e that the spiral 
arms are necessarily a d ensity wave llLin Rr, Shiil If 9641: 
iBinnev fc Tremaiuelll 987^ , and a few even assume a par¬ 
ticular origin of the spiral arm p erturbation. 

For example, iLin et al, ] iTT^ argued that the only 
solution possible for spiral arms are those described by 
the aforementioned density wave theory. Under small 
perturbations, a bes t fit H p = 11 — 13 km sec“^ kpc“^ 
was obtained. lYuanl 11196911 found that the observed 21- 
c m line profile is consistent with the Dp range found 
bv iLin ^ ~al\ lll969ll . However, iNelson fc Matsiida lll977ll 
compared the shock solution for the profile of the gas 
perturbed by a spiral wave to the 21-cm line observed in 
the direction opposite to the center of the galaxy, and 
derived that flp w 20 km sec“^ kpc~ ^. The explanation 
for th e differen t resu l ts, acc ording to INelson fc Matsud"^ 
<|1977fl . is that lYuail lll969ll used a l inear theory for the 
gas motion, and that the ILin et ajl lll969ll results could 
have been modified by the existence of a high velocity 
halo. _ 

iBissantz et al\ ll2003ll modeled the Milky Way as a 
barred galaxy and allowed separate pattern speeds for 
the bar and for the spiral structure. By using hydro- 
dynamic simulations and comparison to the observed 
V — i maps in CO, they showed that the spiral arms 
penetrate the bar’s corotation region, such that within 
the same region of the bar, there are two different 
pattern speeds. The bar’s pattern speed obtained is 
« 59 ± 5 km sec“^ kpc“^, and that of the spiral arms 
is « 19 ± 5 km sec“^ kpc“^ (while assuming that the 
sun is located at 8 kpc from the galactic center). 

A second method which still relies on the assump¬ 
tion that the spiral arms are described by the den¬ 
sity wave theory, is through the study of the kine¬ 
matics of various populations, which orbit the galaxy 
in orbits modified by the spiral wave. Because the 
method relies only on kinematics, the determination of 
the age of the spiral tracers is not important as it is 
for other methods (e.g., when using star c l uster gra¬ 
dients as is done here). iComeron & Torr^ l|199m ob- 
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tained Afip = — Up = 10 ± 5 km sec“^ kpc“^, thus 

rip = 16 ±5 km sec“^ kpc“^, and a pitch angle of i ss 8° 
for a sun located at 8.5 kpc from the galactic center, with 
an angular speed of 220 km sec“^. 

Another method which still relies on the density 
wave th e ory, i s the identification of resonance features. 
laordonl iITotI determined the existence of a disconti¬ 
nuity in the CO abundance as a function of radius, at 
the region of i? = 4 kpc, and identified it with the in¬ 
ner Lindblad resonance. This result is consistent with 
r2p = 11.5±1.5 km sec“^ kpc“^, for a two spiral-armed 
Milky Way. Gordon’s analysis was for the region of 
2 kpc < i? < 10 kpc. 

In a second type of methods, the physical origin of the 
spiral arms is not assumed. Instead, objects which are 
linked to spiral arms and can be “dated”, are analyzed. 
Specifically, either the age gradient or the birth place 
of particular objects can be studied to re construct the 
dynam ics of the spiral arms. For example, IPaloiis et al\ 
fioTl fitted the place of birth of open clusters to a mod¬ 
eled Milky Way based on 21-cm observations. They ob¬ 
tained two possible values: Up « 13.5 km sec“^ kpc“^ 
and rip « 20 km sec“^ kpc“^. 

Interestingly, various authors actually found 
multi-pattern speeds co -existing in the galaxy. 
IRatav fc GottesmanI II2002I1 presented a two dimensional 
model for multiple pattern speeds in barred galaxies. 
Their resul ts are in agreement with Sellwood fc Soark^ 
lll988ll an d IRautiainen fc Said ||199E) for by red galax- 
ies. iKuiikenI II1996I1 . IPicaud fc GarzoiJ II2002I1 and many 
other observations support the possibility of bar in the 
Milky way, which may explain the origin of multiple 
harmonics in the arms. 

To summarize, the current determinations of flp are 
inconclusive. In this paper, we will try to estimate the 
value of rip using the birth place method, for young clus¬ 
ters (age < 5 X 10^ yr), without limiting ourselves to 
assumptions of a particular theory, such as the density 
wave theory. We will show that the origin of the large 
confusion rests in the existence of at least two spiral sets 
besides the Orion armlet, within which the solar system 
currently resides. 

The last parameter characterizing the spiral arms is 
their amplitude. We will not dwell on this point here, 
but a more elaborate analysis aimed at finding the fi¬ 
nite amplitude based on the same set of open clusters, is 
underway. 

2. SPIRAL STRUCTURE ANALYSIS 

Here, we present the methods used to map and study 
the dynamics of the Milky Way spiral arms i n our galac¬ 
tic vi cinity. In our an a lysis, we used the IDias et 'al\ 
ll2f)nl and lLoktin et al\ lll994H databases of open clus¬ 
ters. These include the galactic coordinates and age es¬ 
timates for open clusters at typical distances of up to 
a few kpc from the sun. We used lOlling fc Merrifiel^ 
lll998li for the rotation curves. 

We work in cylindrical galactic-centered coordinates, 
r, (j) and z. The solar system’s position in the Milky 
Way has not been accurately ascertained yet, conserva¬ 
tively, it ranges between 7.1 kpc an d 8.5 kpc, based on 
veloc ity-longitude maps of GO data llOlling fc Merrifield 
[T99|). Neyertheless, recent determinations based on the 
dynamics of stars around the central black hole do giye a 


tighter estimate. For example, lEisenhauer et al\ (1200311 
find a distance of Tq = 7.94 ± 0.42 kpc based on the ob¬ 
served dynamics of stars ar ound the central black hole , 
while using Hipparcos data, iStanek fc Garnavichl (1199811 
find r© = 8.2 ± 0.15stat ± 0.15sys kpc. We will be con¬ 
servative henceforth and consider that the larger range 
of 7.1 — 8.5 kpc is possible, but keep in mind that the 
higher values are more likely. 

Similarly, the reasonable solar rotational velocities Vq 
range between 184 and 200 km sec“^ for a solar system at 
r© = 7.1 kpc, and 220 — 240 km sec“^ for r© = 8.5 kpc. 
We carry out all our analyses while assuming five possible 
locations and velocities. The parameters of these five 
different possible configurations for the sun orbiting the 
galaxy are r© = 7.1 kpc and r© = 8.5 kpc with the 
lower and upper u© of the two aforementioned ranges, 
and r© = 7.8 kpc with u© = 202 km sec“^ as a nominal 
case. 

We conduct our analysis in a frame of reference rotat¬ 
ing with the spiral pattern speed, Hp. In this frame, the 
differential angular velocity of each cluster is given by 

AH = H(r) — Hp. (1) 

Furthermore, we consider only relatively young clus¬ 
ters, having ages satisfying t < 5 x 10^ yr. This implies 
that we preform our analysis only on clusters which did 
not have time to move too far away from the arm^. 

The most detailed v — £ maps of the Milky Way are 
the lDame et al\ (1200111 observations of molecular gas us¬ 
ing CO as a tracer. In this data set, one can clearly 
identify three different arms: The Carina arm which 
is located inwards to the solar galactic radius r©, the 
Perseus arm, externally located relative to r©, and the 
“outer arm” beyond the Perseus one. Our Sun is located 
in the Orion arm, which is p robably only an “armlet” 
(iGeorgelin fc GeorgelinllT976li . 

We assume the most general case in which the ob¬ 
served arms are not necessarily part of the same set of 
spiral arms. As we shall demonstrate, this is necessary 
considering the multiple pattern speeds which appear to 
be present. Therefore, we analyze each arm separately, 
and choose to work only with young star clusters (hav¬ 
ing t < 5 X 10^ yr, as explained above). We perform our 
analysis on the Carina arm, the Orion armlet and the 
Perseus arm. 

We work in the frame of reference of the arms. In this 
frame of reference, the birth location is given by^: 

(l^birth = 0 - [f^(r) - Hp] ■ t, (2) 

wher e t is the age of the clust ers, and H(r) was taken 
from lolling fc Merrifieldl (Il998fl . We find a best fit for 
three parameters: Hp, the inclination angle tan(i) = 

I d(ln r)/c?(/)I and the intersection point with the (/)-axis. 
For each Hp we find the best fit, according to a weight 
function defined below. 

^ Typically for clusters, AH < 13 km sec“^ kpc“^ (see 
therefore, the time period for a rotation in the spiral arm frame of 
reference is > 6 x 10® yr. For a four arm set, clusters will cross 
arms every ^ 1.5 x 10® yr. Hence, in < 5 x 10^ yr, typical clusters 
will not be displaced much from the arm. 

® This neglects the finite amplitude of the spiral density wave. 
In our ongoing research, we alleviate this assumption by including 
the epicyclic motion arising from the finite spiral arm perturbation. 
Also, as young objects, the peculiar motion of the open clusters is 
negligible. 
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2 .1. Statistical Methods 

We continue now with the description of the algorithm 
we employ, and in particular, the statistical methods 
used. We begin by assuming a logarithmic arm profile 
llR.oberts et a/jlT97^ . given by 

—— gf^0birth ( 3 ) 

Here a is related to the pitch angle i through the relation 
i = tan“^(a). In other words, a describes the slope of 
the spiral arm in the Inr vs. (purth plot, while b is the 
intersection point with the In r axis. 

We wish to find the best ht for a logarithmic arm. 
Thus, we define a weight function W with which we eval¬ 
uate the goodness of the fit. The definition we chose is 

^ / A- \ 

>V(a,6) =-^max f [l --^],0j (4a) 

i—1 ^ ^ 

Ai = I luA - {a(j)birth,i + &)|, (4b) 

where N is the number of data (i.e., cluster) points, and d 
is a scaling factor describing the typical width of the arm. 
The argument in the summation is positive for a cluster 
inside the arm, and larger for a cluster closer to the cen¬ 
ter. The condition Ai < d = 0.05 enables us to disregard 
data points which do not lie within ~ 700 pc. In other 
words, we effectively consider the width of the arms to 
be 700 pc. The advantage of such a choice for W over a 
more standard dehnition (such as least squares), is that 
under a standard definition such as the sum of distances 
squared, outliers are very important for the minimiza¬ 
tion. However, in our case, it is most likely that outliers 
are simply unrelated to the arm and should therefore be 
disregarded altogether. This is because the definition is 
such that when In(ri) = a(l)birth,i + b, the argument in 
the summation obtains its maximum value of 1, while 
for |ln(ri) — {a4>birth,i + b)\ > d the argument obtains its 
lowest value 0, such that the cluster is disregarded. 

The weight function W is defined using the distances 
in the Inr axis, since in any usage of i^birth (as defined in 
eq. 121 involves the age of the cluster which has a higher 
error. The best ht for a given Hp is then obtained by 
minimizing W. A best ht is subsequently obtained for 
different Hp as well. 

We considered two minimization methods for W; these 
were the metrop olis simulated annealing method (e.g. 
iPress et n./.ll1 9Sfi|l and a “brute force” method, where we 
loop over the whole parameter space with a hnite res¬ 
olution. The advantages of the former method are the 
smaller CPU requirements and higher accuracy reached. 
The disadvantage, on the other hand, is the uncertainty 
of whether an obtained minimum is indeed the best ht. 
For this reason, we also performed a “brute force” min¬ 
imization by looping over parameter space with a h- 
nite resolution, to ensure that the simulated annealing 
reached the absolute minimum. 

The brute force method also proved useful because the 
additional local minima could be identihed. As we shall 
see below, an additional minimum is related to an addi¬ 
tional set of arms, coexisting in the data. 

Since we have no or very poor knowledge of the sta¬ 
tistical properties of the weight function (eq. 0J, nor of 
the statistical or physical variance in the cluster data. 



Fig. 1. — Weight function vs. the pattern speed, Qp, for the lower 
(Carina arm). The five graphs represent different assumed vq and 
vq. Each graph is obtained by minimizing the weight function 
while fixing the value of flp. The two minima in W demonstrate 
that there are probably two sets of spiral arms. 


we used the Bootstrap Method (e.g., IPress et~a!\\198(t\ to 
estimate the confidence levels around the obtained min¬ 
ima. This method allows us to evaluate the errors of the 
obtained results using the data itself, hence its name. 

Specihcally, we estimate the error in the best fit as 
follows. Our best estimate for the best ht parameters 
is obtained by minimizing the weight function W using 
the original cluster data. To estimate the error on the 
best ht parameters, we then build a set of different real¬ 
izations of the data, which differ from the original set by 
having a random 37% ^ 1/e of the clusters replaced with 
random clusters chosen from the original set. We then 
apply the minimization procedure to each new data set. 
The estimated best ht parameters in each realization are 
expected to vary from the best estimate (based on the 
original set) by ^ Icr. Thus, we can use the ensemble of 
modihed sets to calculate the errors. 

3. THE DYNAMICS OF THE CARINA ARM 

We now proceed to apply the above methods to the 
Carina arm. We begin by calculating the best ht for the 
spiral arms as a function of the assumed Up. We plot the 
minimized weight function W as function of Up and hnd 
a bimodal pattern (see hg. ^). This can be interpreted 
as two different sets of spiral arms in the Milky Way. 

A second apparent feature of hg. ^ is that minimiza¬ 
tion hts having a larger solar-galactocentric radius yield 
overall lower minima. This strongly suggests^ that cq > 
8 kpc. Interestingly, it is consistent with various recent 
determination of the distance to the galactic center, as 
mentioned in m 

The First Set: From the two minima we find (see fig.^, 
the minimum corresponding to the lower pattern speed in 

This assumes that the spiral ar ms are logarithmic s pirals, as 
favored in the density wave theory iRoberts fl/jfl975l) . If this 
assumption is alleviated, lower galactocentric radii could still be 
possible. 
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the angular velocity “spectrum” appears to be a some¬ 
what better fit. The results for the different fits to a 
logaithmic arm, while assuming different solar galacto- 
centric radii and velocities, are depicted in fig. [21 The 
actual values of this best fit (henceforward, the “slower 
Carina arm”) are summarized in table Q 
The Second Set: The Carina-arm cluster data include 
a second statistically significant minimum. The signifi¬ 
cance can be inferred from the fact that all realization 
which have 1/e of the clusters randomly exchanged with 
clusters from the data set still exhibit this minimum, 
albeit with a slightly modified best-fit parameters. We 
therefore interpret the data as two existing spiral arms 
with two different velocities which happen to spatially 
overlap at present. The actual best-fits for this solution 
(henceforward, the “faster Carina arm”) are plotted in 
fig- El and summarized in tabled 

Note that although the pattern speed here is faster 
than in the first set, it is close to the orbital angular 
velocity of the sun, implying that the “synodic” angular 
velocity of the faster set is relatively sma ll. 

As summarized in l[T]and elaborated in iShavi l (l2nnl . 
the determinations of Qp appear to divide into three dif¬ 
ferent ranges of pattern speeds. The two values obtained 
for rip in the Carina arm (see fig.^ nicely agree with two 
of the three different ranges obtained by other authors. 
We will elaborate on the implications of the results in 
the discussion. 

4. THE DYNAMICS OF THE PERSEUS ARM 

When compared with the Carina arm, the Perseus arm 
includes less data points. Consequently the results we 
obtain are not as conclusive as we could hope for, though 
they are still statistically significant. The values found 
are summarize in table 111 and fig. El 

The pattern speed obtained for the Perseus arm is 
somewhat higher than the lower Carina-arm pattern 
speed, though not totally consistent if one considers only 
the statistical error. This could arise from unaccounted 
systematic errors, as will be elaborated in the discus¬ 
sion. Moreover, if we repeat the analysis while fixing 
the Perseus arm pitch angle with the somewhat higher 
value obtained for the slower Carina arm, then the best 
fit obtained for the Perseus arm has a ~ 7r/2 phase shift 
relative to the slower Carina arm. This strongly suggests 
that the two arms are indeed related to each other as part 
of an m = 4 set. If we repeat the analysis while fixing 
rip for the Perseus arm with the value of the slower Ca¬ 
rina arm, we again obtain a fit which has a ~ 7r/2 shift 
relative to the slower Carina arm. 

5. THE DYNAMICS OF THE ORION ARM 

The Sun is located in the vicinity of the Orion arm. 
The Orion arm itself is often regarded as an armlet, 
since it does not appear to be part of the Cy ina and 
the Perseus arms llOeorgelin fc CeorgelinlilOT^ . Never¬ 
theless, it is a site of star formation and has the features 
of a density wave. Irrespective of whether or not it is 
a full hedged spiral arm, or just a “small” pertrubation, 
we can estimate the pattern speed associated with it. 

Looking at fig. El we can identify a high pattern speed 
for the entire range of the assumed solar galactocentric 
radius and angular velocity. We should point out that 
for the 8.5 kpc and the 7.8 kpc configurations, we find 


two different patterns, with very similar pattern speeds 
but somewhat different pitch angles. One may argue 
that this indicates multiplicity is the spiral structure, but 
since the two pattern speeds resemble in value, it is not 
a strong argument. This is because we only have rela¬ 
tively few cluster points to represent the Orion arm, such 
that the very few clusters at large radii can signficantly 
offset the pitch angle determination. Moreover, for a 
large assumed pitch angle for the Orion arm, the Orion 
and Perseus arms coincide, in which case Perseus clusters 
can accidentally be included in the pitch angle determi¬ 
nation. Thus, we can only be certain of the existence of 
one set. The fit having the lower pitch angle (and lower 
pattern speed) is more favorable, while the second solu¬ 
tion could be an artifact, because of the aforementioned 
coincidence. 

Irrespectively, we summarize our results for both fits 
of this arm, but the multiplicity should be taken with a 
grain of salt. An example fit to the Orion arm is por¬ 
trayed in fig. El while the rest of the data is summarized 
in tabled 

6. DISCUSSION 

In our analysis, we found two sets of spiral arms co¬ 
existing in the vicinity of what appears to be the sin¬ 
gle “Carina arm”. In the “velocity spectrum” diagram, 
the lower rip fit is more prominent, with lower val¬ 
ues of the weight function W (see figure EJ. It corre¬ 
sponds to a pattern speed of flp, Carina, i = ± 

l.lstat km sec“^ kpc“^, and a pitch angle of *carina.i = 
15.8lQ gj,j^^ ± The systematic error only includes 

the uncertainty arising from the inaccurate determina¬ 
tion of the solar galactocentric radius and velocity. Other 
systematic errors could certainly exist, though hard or 
impossible to properly account. For example, the clus¬ 
ters we use are biased towards lower distances, implying 
that a higher weight is given to clusters on one side of the 
arm. In the case of the Carina arm, it is the older clusters 
which are favored. In the Perseus arm, it is the oppo¬ 
site. A second systematic error arises from our implicit 
assumption in eq. (21 that the amplitude of the density 
wave is negligible. Alleviation of this assumption is done 
in a paper in progress. 

Since the measurement is essentially based on the 
speed of the star clusters relative to the orbital frequency 
at our galactic radius, a more accurate result is obtained 
if we fit for An = IIq — flp, such that the systematic 
error in Hq is less important. We find: 

Hq Up, Carina,! — I0.6_q i l.Iaiai km SeC kpC 

For comparison, the relative pattern speed Afl can be 
measured using the periodic variations in the cosmic ray 
flux history as recorded in Iron meteorites. Th is gives 
An = 10.5 ± 1.5sys ± O.Sstat km sec“^ kpc“^ llShavivI 
l2finl . Alternatively, it can be compared to the period¬ 
icity obtained through the climatic effects of the vari¬ 
able cosmic ray flux. This geological determination re- 
sults with Ar i = 10.4 ± 1.5sys ± 0.35stat km sec“^ kpc“^ 
(jSh avi 1 12001 . Thus, the lower pattern we find for the 
Carina arm is in agreement, within the measurement er¬ 
ror, with both a direct and indirect determination of 
the cosmic ray flux variability. Note that the very good 
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r0. The solid line depicts the best fitting logarithmic arm, while the dotted lines denote ±700 pc strips delineating the clusters which are 
considered as part of the spiral arm. 



agreement between the results obtained here and the me- 
teoritic or geological determination should be considered 
a coincidence. The agreement is less impressive if one 
consolidates the results of the slower Carina arm with 
the Perseus arm, as discussed below. 

The rip, Carina ,1 pattern speed places the sun in the 
vicinity of the inner (1:4) Lindblad reson ance (se e hEi.Ell . 
This is consistent with the analysis of iQuilli^ ll200,‘lll . 
who modeled stars in the solar vicinity and found that 
the observed dynamical signatures of chaos can be ex¬ 
plained if the solar system is in the vicinity of an in¬ 
ner Lindblad resonance (either the 1:4 resonance for a 4 
armed set or the 1:2 resonance for a 2-armed set). 

With reasonable certainty, we can also state that a 
weaker arm, possibly part of a larger set, is overlapping 


the previously described arm at the same radial range, 
but with a higher pattern speed. The second pattern 
speed we obtained for the Carina arm is Dpcarina ,2 = 
29.8l2^2sys =*= 2.1stat km sec“^ kpc“^, with a pitch an¬ 
gle JCarina .2 = 21.ltJ ± 1.3°t„j. The Statistical error 
here is smaller than for the faster arm, with the reason 
being the smaller dispersion associated with the smaller 
relative speed. The relative pattern itself is 

Hq f2p,Carina,2 — 2.T ib Igys i O.Qstat km SeC kpC 

( 6 ) 

Again, we find that the relative pattern speed is better 
constrained than the absolute value. 

The value of np,Carina ,2 places the observed outer ex¬ 
tent of the Milky way’s spiral structure in the vicinity 
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TABLE 1 

Spiral pattern speeds for given vq and vq 


r© [kpcj 

7.1 

7.1 

8.5 

8.5 

7.8 

VQ [km sec'll 

184 

200 

220 

240 

202 

The First Carina Arm 

Up [km sec ^ kpc ^1 

15.8 ±0.9 

17.9 ±0.7 

15.6 ± 0.9 

16.9 ± 1.8 

15.3 ± 1.6 

i 

16.6 ± 1 . 2 “ 

16.2 ± 1.3“ 

15.6 ±0.4“ 

14.9 ± 1.9“ 

15.7 ± 1.5“ 

The Second Carina Arm 

Up [km sec ^ kpc ^1 

29.7 ±0.8 

31.8 ±0.6 

27.6 ± 0.3 

30 ±0.5 

28.4 ±0.6 

i 

20.5 ± 1.9“ 

21.1 ± 1.7“ 

21.4 ±0.9“ 

21.5 ± 1“ 

22.4 ± 1.1“ 

'The Perseus Arm 

Up [km sec ^ kpc ^1 

20.0 ± 1.7 

21.9 ±1.3 

18.9 ± 1.9 

21.0 ± 1.8 

18.9 ± 1.3 

i 

10.2 ± 1.5“ 

9.4 ±0.4“ 

8.5 ± 0.5“ 

8.5 ±0.5“ 

8.5 ±0.5“ 

'The Orion Arm 

Up [km sec ^ kpc ^1 

27.9 ±0.5 

30 ±0.5 

27.7 ±0.5 

30.2 ± 1.6 

27.4 ±0.5 




31.5 ± 0.5 

33.4 ±0.5 

31.8 ±0.5 

i 

10.6 ± 0 . 6 “ 

10.5 ±0.6“ 

8.5 ± 0.7“ 

8.5 ±0.5“ 

9.2 ±0.4“ 




21.3 ± 0.9 

21.3 ±0.9“ 

21.7 ±0.7° 



Fig. 4.— The circular frequency H (solid line) and the 1:2 (dash- 
dotted line) and 1:4 (dotted line) Lindblad resonance frequencies 
as a function of the galactic r adius. Given tha t spiral arms are 
observed out to router ~ 2 rQ iBlitz et fl/Jll983l . this is probably 
the outer Lindblad resonance. Nevertheless, the outer Lindblad 
resonance could in principle also be at larger radii. The four arms 
could be part of one set with flpi. In this case, the outer extent of 
the spiral arms, should be identified not with the outer Lindblad 
resonance, but with the co-rotation. Interestingly, here the inner 
Lindblad resonance is near tq. We also see that as most two of 
the arms could be part of an m = 2 set with In this case, 

the outer extent of the spiral arms is just near the expect outer 
Lindblad resonance. 


of the outer (1:2) Lindblad resonance. This also implies 
that we are just outside co-rotation for this spiral set. 

Our results are in agreement with the two ranges 
of values often found in the literature, of ~ 
16 — 20 km se c~^ kpc ~ ^ and ~ 30 km seckpc“^. 
For instance, llvano 3 <1198311 found flp = 16 — 

20 km sec“^ kpc“^, using cluster age gradients. His 
result agrees with our first set, as w ell as with 
iNelson fc Matsud"^ <1197111 . lEfremovI <1198311 . and others. 
The second number is in nice agreement with previous 


determ inations as well. For example, iFernandez et all 
<l20(rill found O p = 30 ± 7 km sec ^ kpc ^, while 
IMishurov fc Zeninal <1199911 placed the solar system at 
AR > 0.1 kpc, relative to the co-rotation radius. 

From linear density wave theory, spiral arms are ex¬ 
pected to exist only between the Lindblad resonances 
depending on the particular pattern speed. Observation- 
ally, a clear 4-armed spiral set appears to exist at least 
from our galactocentri c radius out to about twice 
this radius. Tout ~ 2 r 0 <IBlitz et ol jiMa. 

When combining the above theoretical argument with 
observations, the outer extent of the faster set should end 
significantly inwards of the actual observed outer radius 
Tout, if the fast arm is part of a 4-armed set. On the 
other hand, it could extend out to Tout , if the arm is part 
of a 2 -armed set (see fig. El). In other words, the 4-arms 
observed out to rout cannot all belong to the faster set. 
Thus, one possibility is that both the slower and faster 
sets are 4 armed, in which case the slower set extends out 
to Tout while the faster set necessarily ends significantly 
inwards of rout- A second possibility is that the slower set 
is still 4-armed but the faster set is 2-armed. Here, the 
observed 4 arms are still part of the slower set, however, 
superimposed on them, we have a faster 2 -armed set that 
could extend as far out as rout , but could also end further 
in. 

The inner extent of the observed spiral arms can be 
used to constrain the patten speeds from below. Inspec¬ 
tion of fig. El reveals that the inner extent of the slower 
set should be roughly at our galactic radius. In fact, the 
nominal lower limit is just outside the solar circle. Since 
we know the sun does pass through the arms of this set 
(e.g., from meteoritic or geological data lShavivll2003ll . a 
very small discrepancy of less than a few 100 pc exists. 
This should no be a source of major concern for several 
reasons. 

First, our knowledge of Qp which determines the inner 
Lindblad radius suffers from both systematic and statis¬ 
tical errors. Increasing by its statistical error, for 
example, shifts the inner extent by more than 0.5 kpc. 

Other interesting points to consider is that the sun per¬ 
forms epicyclic motion around a center located roughly ~ 
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Fig. 5.— A plot of Inr vs. the birth angle </ibirth of clusters and the best fit for the Perseus arm, as described in fig.E] 



Fig. 6 .— The weight function W vs. the pattern speed flp for 
the Orion arm. The five different graphs represent the different 
assumed location and velocity of the sun with respect to the center 
of the galaxy. Either one or two nearby minima exist, for flp 
^ 0 . Namely, the solar system is near the co-rotation radius of 
the Orion arm. In the inset we plot W vs. Qq — Qp. Evidently, 
the method is more accurate in determining the relative pattern 
speed of the arms. Similarly to the Carina arm, we again find that 
configurations with tq '^8 kpc give consistently better fits. 


0.5 kp c outside the current solar circle le.g. lMatese et ali 
allowing the sun to pass through the spiral arms 
even if it is currently just at their inner edge. 

Near the resonance, the response of both the gas 
and stars to the spiral wave perturbation is extreme. 
This is expected to give rise to interesting nonlin¬ 
ear effects. Moreover, in several numerical simulations 
and observation, spiral arms have been observed to 
exceed the in ner, outer or both Lind blad resonance. 
For example, iSemoere fc R.ozaj lll997^ have Ha de- 



the Orion arm, while considering one configuration for the galaxy 
(similar to fig. 1 ^. 



Fig. 8 .— A histogram showing Q,p of the best fits for the Orion 
arm, revealing the notable bi-modality. 
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projected images of the galaxy NGC 157, demonstrat¬ 
ing the existence of spiral arms beyond the pr e dicted 
oute r Lindblad r esona nce. iR.egan fc Mulcha^ lll999ll 
and iLaine et a,l\ <ll999tl reported the existence of spiral 
patterns which exist inside the inner Lindblad resonance 
and connected to the outer pattern. This prompted 
lEnglmaier fc ShlosmanI (1200(11 to analyze the formation 
o f spiral pa t terns inside the Lindblad resonance. 

ISellwoodI (1200011 reached similar conclusions about spi¬ 
ral arms beyond the Lindblad resonances, but based 
his analysis on numerical simula tions. Furthermore, 
iGoldreich fc Tremaind (|1978l[T97^ studied the propaga¬ 
tion of waves in spiral galaxies. They argued that short 
density waves are not restricted to any region (whereas 
long density waves are restricted to exist between the 
Linndblad resonances). This is especially the case, if the 
origin is a rotating bar. 

To summarize the “Carina arm”, we can conclude that 
it is a superposition of two arms with two different pat¬ 
tern speeds. The slower arm is definitely 4-armed in 
structure and dominates the outer parts of the galaxy. 
On the other hand, nothing definitive can be concluded 
on the number of arms in the faster set, which dominates 
the galaxy within the solar circle. 

We can also conclude, more generally, that there is a 
real physical reason for the previous discrepancies be¬ 
tween the different pattern speed determination—any 
pattern speed analysis which does not allow for the pos¬ 
sibility of more than one set should be considered cau¬ 
tiously. 

For the Perseus arm clusters, we find only one 
clear pattern speed having flp^perseus = 20 . 0 j'ji[' 2 sj,s ^ 
l.Gstat km sec“^ kpc“^, with a pitch angle of iperseus = 
9-Ot n This result is consistent, for example, 

with iMartos e.t nJ\ l|2nfl4f) . though a more careful compar¬ 
ison reveals some small discrepancy as described above. 
This could arise from yet unaccounted systematic errors. 

If we consider that the observed clusters in the Perseus 


arm span the range between ^ to ~ l.Sr©, and that 
there is no sign for the r 2 p ,2 pattern in the data, we 
can further extend our conclusion on the nature of the 
spiral arms, which we obtained using the Carina arm. 
Specifically, if the faster pattern is 4-armed, its outer 
extent should be within corotation (and not, for example, 
the outer 4:1 Lindblad resonance), since corotation with 
r2p^2 takes place in the vicinity of tq. If the faster set is a 
2 -armed pattern, we do not expect an arm superimposed 
over the Perseus arm, and no conclusion can be reached 
regarding its outer extent. 

Concentrating now on the Orion arm, we obtained that 
the best fitting arm has a pattern speed of Op^orion,i = 
28.9lj;2sy5 ± O.Sstat km sec“^ kpc“^, with a pitch an¬ 
gle of iorion.i = ± 0.6°t„(. Again, the system¬ 

atic error arises from our limited knowledge of the solar 
galactic radius and velocity. If we look at AO instead, 
then the systematic error is notably smaller. We find: 
AOorion.l = l-StoJgys =*= 0.7stat km sec ^ kpc ^. 

However, under several assumption about the so¬ 
lar location and velocity, we find a second pos¬ 
sibility, consistent with AOorion.i = 5.6to;4sys T 
O.Sstat km sec“^ kpc“^, and a pitch angle of iorion.i = 

21.4t°;L,s ± 0.8°taf 

As previously explained, it is not clear whether the 
two solutions indeed imply that two pattern speeds exist. 
One possibility which seems reasonable, is that the sec¬ 
ond AO, with its large pitch angle, is obtained as an ar¬ 
tificial fit where an “arm” extends from the actual Orion 
arm clusters to the Perseus arm clusters, and is therefore 
an artifact. 

Interestingly, the Orion and faster Carina sets have 
similar pattern speeds. It may simply be a coincidence, 
since there is no evidence indicating that they are part 
of the same set, however, there could be some kind of 
dynamic coupling. This is an interesting question to ad¬ 
dress in future research. 
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